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Abstract 

We discuss the unstable character of the solu- 
tions of the Lorentz-Dirac equation and stress 
the need of methods like order reduction to 
derive a physically acceptable equation of 
motion. The discussion is illustrated with the 
paradigmatic example of the non-relativistic 
harmonic oscillator with radiation reaction. 
We also illustrate removal of the noncasual 
pre-acceleration with the introduction of a 
small correction in the Lorentz-Dirac equa- 
tion. 
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Lorentz-Dirac (LD) equation £Q, 



The classical equation of motion of a point 
charge including radiation reaction is the 
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where a dot indicates derivative with respect 
to the proper time r; /** is the external 4-force 
per unit mass and b = J^-. Although the 
LD equation is based in solid physical ground 
(Special Relativity and Electrodynamics) it 
has some unusual mathematical properties 
that put some doubts about the suitability of 
this equation to represent the motion equa- 
tion of a charged particle with radiation reac- 
tion. Since the LD equation is a third order 
ordinary differential equation, not only the 
initial position and velocity, but also the ini- 
tial acceleration must be given to specify a 
unique solution. Moreover, one has to de- 
termine with an infinity accuracy the initial 
value of the acceleration that eliminates the 
part of the solutions that grow exponentially 
in time ("runaway" solution). 
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Recently, Spohn [2] found that the physical 
solutions of the LD equation lie on a critical 
manifold in phase space that consists exclu- 
sively of repelling fixed points. On this criti- 
cal surface the motion is free of runaway so- 
lutions and is governed by an effective second 
order equation. If an initial condition slightly 
off the critical surface is given, the solution 
grows exponentially fast. This lack of stabil- 
ity of solutions with respect to small varia- 
tions of initial conditions means that the LD 
equation does not represent a mathematically 
well posed problem [3] and in consequence its 
applicability is limited. 



Difficulties also arise in the numerical so- 
lution of the LD equation [4 . Even if one 
knows the correct initial value of the accel- 
eration, in the usual forward numerical in- 
tegration the runaway contributions to the 
solution grow extremely fast due to numeri- 
cal errors. As an alternative, Aguirregabiria 
[Sj proposes a method of order reduction in 
which Eq.(JTJ) is substituted by a second or- 
der equation of motion with no runaway so- 
lutions. The explicit form of this second order 
equation cannot be obtained in general, so a 
series of successive approximations x M = 



(n = OA, . . .) is constructed: 
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We note that in this approach only the usual 
initial conditions (position and velocity) are 
required and the equations for a finite n are 
stable with respect to variations of initial con- 
ditions. 

Rohrlich , inspired on Landau [7] , argues 
that the original LD equation should be re- 
placed by Eq. (j3J) as the exact equation of mo- 
tion of a radiating point particle, since it is 
of second order, satisfies the principle of iner- 
tia, and the radiation reaction term vanishes 
in the absence of an external force. This pre- 
scription is based on physical arguments, but 
only the confrontation with experimental re- 
sults will establish its validity. For proposals 
of experimental tests, see [2] . 

The purpose of this letter is to illustrate 
the above discussion with the significant ex- 
ample of the non-relativistic one-dimensional 
harmonic oscillator with radiation reaction. 
This problem has an analytical solution that 
can be used to construct the effective sec- 
ond order equation of motion on the above 
mentioned critical surface. The successive 
approximations scheme (0)-© and solutions 
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from direct numerical integration of the non- 
relativistic LD equation can be confronted 
with the exact solution. The study of this 
example will show us the limitations of the 
LD equation 

to describe the motion of a charge particle 
interacting with its own radiation (back reac- 
tion). Also we shall discuss the addition of a 
small correction term that eliminates the pre- 
acceleration in the solutions to the LD equa- 
tion [8 . In particular we discuss, in some 
detail, the case of a constant force that acts 
on the charge particle during a finite time. 

The one-dimensional non-relativistic mo- 
tion of an electron subjected to a harmonic 
force is governed by the Abraham-Lorentz 
limit of the LD equation, 

d cc o - d cc , 

whose exact solution is given by 

x(t) = de ait +e a2rt (c 2 cos(a 2 it) + c 3 sin(a 2 ^)) 

(6) 

The constants ci, c 2 and c 3 depend on the 
initial conditions and 



^From Eq.fjHJ) we find the expressions for 
the velocity and acceleration: 

x(t) = a.\C\e ait + (c 2 a 2r + c 3 a 2i )e a2rt cos(a 2i t) 
+ (c 3 a 2r - c 2 a 2i )e a2rt sm(a 2i t), (12) 
x(t) = a\cie ait + [c 2 (a 2r - a 2j ) + 2c 3 a 2r a 2 i] x 
e a2rt cos(a 2i t) + [c 3 (a 2r - a 2 2i ) 

-2c 2 a 2r a 2i ] e a2rt sm(a 2i t). (13) 

In terms of the initial conditions x(0) = xq, 
x(0) = xq and x(0) = xq, the constants ci, c 2 
and c 3 are given by 



ai = ^ (1 + P 1/3 + /T 1/3 ) , (7) 



a 2r 



3b 
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^ (f 3 - r 1/3 ) 



(9) 
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8 + 108cj 2 6 2 + 12v / 3cj6v / 4 + 27cj 2 6 2 . 



ci 



x (a 2r + afj + x - 2a 2r x 
(at — a 2r ) + a 2i 



(14) 
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x «i (oil - 2a 2r ) - x + 2a 2r x 

C2 = ^—— 2 > ( 15 ) 

[ai - a 2r ) + a 2i 

c 3 = — [xqoci {a\a 2r — c? 2r + a 2i ) + 
+x (ai — a 2r ) + xq (a 2r — a\ — a 2i )] / 

a 2i [{a x - a 2r ) 2 + a 2i ] . (16) 
Using Eq.© and Eq.flEJl, Eq.(JT3J) can be 

Ccist dbS 

x{t) = ci [(oti - a 2r ) 2 + aji] e ait + 2a 2r x 

-(a 2 2r + a 2 2i )x. (17) 

Since a.\ > 0, the term e ait has a runaway 
character, so we must set c\ = to eliminate 
the unphysical solutions. Then Eq.(JT7J) takes 
the form 

x(t) = 2a 2r x — (a 2r + a\i)x. (18) 

This is the equation of motion of a damped 
harmonic oscillator and represents the effec- 
tive second order equation on the critical sur- 
face. In this particular problem the effective 
equation can be obtained exactly. 
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Figure 1: Acceleration as function of time 
with bu = 2tt/7 and initial conditions x(0) = 
1, x(0) = (solid curve), x(0) = ±1 x 10" 10 c 
(dashed curves) and a;(0) given by Eq. ffTHj) . 
The curves were ploted using the analyti- 
cal solutions of Eq.(j3J) (1 unit of length = 
2.818 xl(T 15 m, 1 unit of time = 9.399 x KT 24 
s). 



Figure 2: Acceleration as function of time 
with bu = 2n/7 and initial conditions x(0) = 
1 unit of length, x(0) = (solid curve), 
x(0) = ±1 x 10~ 10 c (dashed curves) and 
x(0) given by Eq. (jl8|) . The curves were 
drawn numerically integrating Eq.(jSJ) (1 unit 
of length = 2.818 x 1CT 15 m, 1 unit of time 
= 9.399 x 1(T 24 s). 



The scheme of successive approximations 
applied to Eq.© results in equations similar 
to Eq.JTH) [Sj 

x = — 7„x — ijJ 2 n x , n = 0, 1, ... (19) 

where the coefficients satisfy the recurrence 
relations 

7o = 0, 

^n+l = - bulin, (20) 
7n+i = b(u 2 n - 7 2 ) . 

When bui ^ 0.95 these successive approxima- 
tions will always converge. 

Figure 1 shows the curves of acceleration 
as function of time with x(0) = 1, bu = 2ti/7 



(our units are: 1 unit of length = 2.818 x 
10" 15 m, 1 unit of time = 9.399 x 10~ 24 
s). The continuous curve was obtained from 
Eq.flTHJ) with initial velocity x(0) = 0, the 
dashed curves with the same initial values 
for position and acceleration, but with x(0) = 
±1 x 10~ 10 c and Eq.(JT3J) -(fTHJl. This is equiv- 
alent to take initial conditions slightly off the 
critical surface. The acceleration diverges af- 
ter a few units of time, indicating the insta- 
bility of the critical surface. 

Figure 2 shows the curves of acceleration 
as function of time with bu = 27r/7 and the 
same initial conditions as in figure 1, but 
now Eq.fjSJ) was solved numerically (fourth- 
fifth order Runge-Kutta with adaptive step 
algorithm) . Even though the value of the ini- 
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Figure 3: Acceleration as function of time 
with bu = 2tt/7 and initial condition x(0) = 1 
unit of length and x(0) = 0. The curves were 
plot using successive approximations Eq. (flT?|) 
for n = 20 and n = 30 (1 unit of length = 
2.818 x 10" 15 m, 1 unit of time = 9.399 x 10~ 24 
s). 

tial acceleration given by Eq. (|18|) is used with 
16 significant figures, the runaway contribu- 
tions make also the numerical solution useless 
after a few time units. 

Finally, in Figure 3 the solutions calculated 
with the successive approximations (fT9"j) are 
compared with the exact solution using initial 
conditions x(0) = 1 and x(0) = 0. Although 
the convergence is slow, no instability is seen. 

This example clarifies the uselessness of 
Eq.(^) in practice. Even though, for very 
special cases, some exact solutions of the LD 
equation are known jUj for actual situations 
one must rely on numerical solutions. But, 
to eliminate the runaway components one 
would have to know the exact initial condi- 
tions and to perform the computations with 
infinite precision. The knowledge of the ex- 



act initial conditions imposes a physically im- 
possible condition to be achieved. Also, do 
numerics with infinite precision is not possi- 
ble. Recently, Chicone et al [TU] pointed out 
that high order equations like Eq. (0) must be 
viewed as an intermediate step in the deriva- 
tion of the physically correct, second order 
equation. A possible method to find this 
equations is the use of successive approxima- 
tions. However, there is no guarantee of con- 
vergence. 

We would like to comment another unusual 
feature of the LD equation, the phenomenon 
of pre-acceleration. When the external force 
is an explicit function of time, the solution of 
the one-dimensional non-relativistic version 
of Eq.Q can be written as ^9] 

Ht) = lJ e-^~^ b f(t')dt'. (21) 

It follows from Eq. (j21j) that a nonzero ac- 
celeration exists before the external force is 
applied at t = 0, thus violating cusality. 
In a nice monograph on the LD equation, 
Yaghjian [S] carefully rederives the equation 
of motion for the extended model of the elec- 
tron clS 8b charged insulating sphere of ra- 
dius a. He shows that the multiplication of 
the electromagnetic self force by a correction 
function r)(r) eliminates the pre-acceleration 
from the solution to the original LD equation 
without introducing false discontinuities in 
velocity across r = or spurious delta func- 
tions and their derivatives at r = 0. The cor- 
rection function increases monotonically from 
zero to one in the time it takes light to travel 
across the electron, and it approaches zero 
like t 2 or faster as r approaches zero from 
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the positive (right) side. This modification is 
needed to ensure the validity of the Taylor ex- 
pansion that is used in the evaluation of the 
electromagnetic self force about present time 
(t) of position, velocity and acceleration of 
each element of charge at retarded time dur- 
ing the interval < r < 2a/ c. With this 
correction, Eq. (|2T|l is modified to 



can be evaluated exactly We get 



x(t) 



< t < t 







x(t) 



g b Jt n(t") 



, t < 
dt' ,t> 



(22) 



g 6 ^ t — £q J 

< to < t < t + 2a 
k[l -e-C*!-*)/^ 

t + 2a < t < ti 




b 



A simple example helps to clarify the dif- 
ferences between Eq. (J2T|) and Eq. (j22"|) . Let 

f{t) be a constant force that acts during a 
finite period of time, 

'0 , 0<t<t 

/(/)= {k ./,,-. /• /; (23) 

,t>t 1 . 
The solution Eq. (|2T|) is given by jH] 

{ ke t/b( e -t /b _ ^ o < t < t 

k[l - e-^-^/ b ] , t < t < ti 

,t>h. 

(24) 

If we assume following form for the correction 
function T)(t): 



T](t) 





(t-tp) 



, < t < t 
j^- ,t <t<t + 2a (25) 
, t > t + 2a, 



where we suppose that t + 2a < t\, Eq.(j22J 



, t > ti. 
(26) 

0- Eq . (|26p shows there is no pre-acceleration 
in the interval < t < to before the exter- 
nal force begins to act. In the limit t — > , 
solution reduces to x(t ) = k, i. e., the 
acceleration equals the external force per unit 
mass when the external force is first applied, 
a result that Yaghjian shows is valid in gen- 
eral. Figure 4 displays the curves of accelera- 
tions calculated with Eq. and Eq. (}2T)]) for 
to = 2, ti — 5, a — 1, and k = 0.1. When 
r](t) reduces to 1, after t — t + 2a, both so- 
lutions agree. 

To summarize, the LD equation suffers 
from two major problems: pre-acceleration 
and runaway solutions. To have a physically 
acceptable equation of motion two 

approaches have been proposed. In the 
first, we introduce a correction function mul- 
tiplying the electromagnetic self-force term 
in the original LD equation to eliminate the 
pre-acceleration jH] and impose a boundary 
condition to eliminate the runaway solution. 
And in the second by an appropriate proce- 
dure we substitute the third order differential 
LD equation by an effective second order dif- 
ferential equation [7J. Rohrlich advocate 
to take this effective second order equation 
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without correction 
with correction - 



Figure 4: Accelerations calculated with 
Eq.(j2U) and Eq.flHD for t = 2, t x = 5, a = 
l,k = 0.1, (1 unit of length = 2.818 x 1(T 15 
m, 1 unit of time = 9.399 x 10" 24 s). Note the 
elimination of pre-acceleration when the cor- 
rection function is included in the Abraham- 
Lorentz equation of motion. 

as the right equation for the charged particle 
motion. 

Acknowledgments. D.V. thanks capes and 
P.S.L. thanks FAPESP and CNPQ for financial 
support. 



[4] W.E. Baylis, J. Huschilt, Phys. Rev. D 13 
(1976) 3262. 

[5] J.M. Aguirregabiria, J. Phys. A: Math. 
Gen. 30 (1997) 2391. 

[6] F. Rohrlich, Phys. Lett. A 283 (2001) 276. 

[7] L.D. Landau, E.M. Lifshitz, The Classi- 
cal Theory of Fields, Butterworth - Heine- 
mann, Oxford, 1975. 

[8] A. D. Yaghjian, Relativistic Dynamics of 
a Charged Sphere, Springer- Verlag, Hei- 
delberg, 1992. 

[9] G.N. Plass, Rev. Mod. Phys. 33 (1961) 
37. 

[10] C. Chicone, S.M. Kopeikin, B. Mash- 
hoon, D.G. Retzloff, Phys. Lett. A 285 
(2001) 17. 



References 

[1] P.A.M. Dirac, Proc. R. Soc. (London) 
A167 (1938) 148. 

[2] H. Spohn, Europhys. Lett. 50 (2000) 287. 

[3] See, for example, E. Zauderer, Partial Dif- 
ferential Equations of Applied Mathemat- 
ics, 2nd. ed., John Wiley & Sons, 1989, p. 
136. 



7 



